The Lorenz number in CeCoIn 5 inferred from the thermal and charge Hall currents 
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The thermal Hall conductivity K xy and Hall conductivity a xy in CeCoIns are used to determine 
the Lorenz number Ch at low temperature T. This enables the separation of the observed thermal 
conductivity into its electronic and non-electronic parts. We uncover evidence for a charge-neutral, 
field- dependent thermal conductivity, which we identify with spin excitations. At low T, these 
excitations dominate the scattering of charge carriers. We show that suppression of the spin excita- 
tions in high fields leads to a steep enhancement of the electron mean-free-path, which leads to an 
interesting scaling relation between the magnetoresistance, thermal conductivity and a xy . 

PACS numbers: 71.27. +a,72.15.Eb,72.20.My,74.70.Tx 



I. INTRODUCTION 

The heavy-electron system CeCoIn 5 exhibits a host of 
unusual electronic properties of current interest. In the 
superconducting state, strong evidence for <i-wave pair- 
ing symmetry has been reported [J 0, H, H, Q. The 
FFLO state involving pairing with unequal spin popu- 
lations in an in-plane magnetic field (H _L c) has been 
proposed Q. The phase diagram in a perpendicular 
field (H||c) has also received wide attention @, Q. In an 
extended region in T-H plane surrounding the supercon- 
ducting (SC) state (labelled I in Fig. Q}i), the resistivity p 
and heat capacity exhibit distinctive "non- Fermi liquid" 
characteristics: p ~ T 0], while the Sommerfeld param- 
eter j(T) ~ logT Q. When H exceeds the boundary 
H S (T), p recovers the Fermi-liquid form p = p + AT 2 
and the unconventional features of j(T) are suppressed. 
The high-field region (labelled II in Fig. QJi) is called the 
Fermi-liquid region. The boundary field H S (T) and the 
upper critical field H C 2 terminate at a quantum critical 
point (QCP) as T — > (the field scale Hz is discussed 
below) . Several parameters characterizing resistivity dis- 
play divergent behavior as the QCP is approached hi. In 
addition, a large Nernst signal is observed in I P.[lfj|. 

To clarify the electronic state in the region I, we have 
measured extensively the in-plane thermal conductivity 
with H||c||z, and the thermal Hall conductivity 
K xy (the Righi-Leduc effect) in crystals with very long 
electron mean-free-path I. In addition, we measured the 
electrical conductivity a = a xx and Hall conductivity 
a xy . The crystals, grown from metallic flux are plate- 
like with the c axis normal to the broadest faces and the 
a axis along one edge (for structure, see Ref. [H[). The 
zero-field k displays a prominent peak below T c (solid cir- 
cles in Fig. [TJa), which arises from the steep increase in £ 
of Bogolyubov excitations in the superconducting state. 
Using the Hall conductivities n xy and a xy , we demon- 
strate the validity of the Wiedemann-Franz (WF) law, 
and then use the Lorenz number to separate the total k 
into its electronic and non-electronic components K e and 
Kb, respectively. From the strong field dependence ob- 
served in Kb, we infer that spin excitations provide the 
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FIG. 1: (a) The phase diagram of CeCoIns adapted from 
Refs. 7, 8], showing some features inferred from our exper- 
iment. The dashed line is the boundary H S (T) between the 
Fermi liquid (II) and non-Fermi Liquid (I) regions reported 
in Refs 0,0- Solid triangles are our values for H S (T). Open 
circles indicate the field scale Hz (this work) above which the 
current ratio lZ a = Z (see Eq. [3]). (b) The T dependence 
of n in zero field (solid circles) showing a large qp peak be- 
low T c . The background term Kb measured at 6 T is shown 
as solid triangles. Below ~8 K, Kb is largely comprised of a 
term k s that is very field dependent (and identified with spin 
excitations) . 



dominant scattering mechanism for the charge carriers 
in the region I. The application of an intense field leads 
to suppression of this scattering channel and a sharp in- 
crease in £. This insight sheds light on the large magne- 
toresistance and the unusual features of the Hall effect. 
We discuss the implications for Cooper pairing in the SC 
region. 
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II. THERMAL AND CHARGE CONDUCTIVITY 
TENSORS 



The thermal resistivity tensor Wij is measured by ap- 
plying a weak gradient (ST ~ 10 mK along length of the 
crystal at 0.5 K). Below 2 K, the strong variation of £ 
with T and H is potentially the largest source of error in 
comparing (measured in finite ST) with (ST = 0). 
We minimized the uncertainties by extensive calibration 
of the RuOz thermometers (glued to the crystal), and 
using very slow field scans (0.1-0.2 T/min. at 0.5 K). 

We emphasize that, because £ attains very large val- 
ues below 10 K, it is necessary to use the full matrix 
inversion to reliably convert the measured tensors Wij 
and pij into their reciprocal conductivity tensors, e.g. 
k-xx — W xx /(Wxx + W£ y ). Experimentally, this means 
that n and a in strong fields cannot be obtained without 
measuring simultaneously the diagonal and off-diagonal 
elements of Wij and p^ (leaving out the Hall tensor ele- 
ments leads to errors in k and a of 30% or more at low 
T). 

Figure [5] compares the field dependences of n (Panel 
a) and the in-plane resistivity p (Panel b) at temper- 
atures from 5 K to 0.5 K with H||c. Above 5 K, an 
increasing field decreases slightly the observed n. With 
decreasing T, however, this trend changes. At 2 K, k 
rises gradually with H. At even lower T (0.5-1.5 K), 
this rising trend becomes firmly established in the nor- 
mal state when H exceeds H C 2 (step in k). In the SC 
region below H C 2, a prominent feature in k is the sharp, 
narrow spike which represents the rapid field suppression 
of the zero-field peak caused by scattering of Bogolyubov 
excitations from vortices [HI, [lH, [3. The spike m k vs. 
H below T c is much larger than previously reported ||. 
This reflects a much longer £ in the present samples. 

The complicated behavior of k(T, H) arises because it 
is the sum of the electronic term n e and a "background" 
term n b carried by charge-neutral excitations (spin exci- 
tations and phonons), viz. 



k(T, H) = K e (T, H) + n b (T, H). 



(1) 



Our main finding is that, in CeCoIns, the charge- 
neutral term Kb(T,H) displays an unexpectedly strong 
H dependence. Its T-profile at H = 6 T is shown as 
solid triangles in Fig. [TJd. 

As previously found @, HH , CeCoIn 5 exhibits a large 
magnetoresistance (MR) (Fig. Oj). The initial positive 
MR (H <3 T) is caused by suppression of superconduct- 
ing fluctuations which we discuss elsewhere [1J|- Our fo- 
cus is on the negative MR that prevails for H >4 T at all 
T below ^30 K. As T decreases towards T c , the negative 
MR becomes pronounced. At 2 K, p decreases by ~2.5 
when H reaches 14 T. Both the sign and magnitude pre- 
clude classical MR associated with the Lorentz force. As 
shown below, the MR results from a steep enhancement 
of £ with increasing field. 
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FIG. 2: (a) Curves of k xx vs. H||c at selected T. The steep 
suppression of the zero-field anomaly produces a sharp spike 
at each T < T c . (b) The magnetoresistance (MR) p vs. H||c 
showing the strong negative MR above 4 or 5 T. 



III. LORENZ NUMBER FROM SCALING OF 

fcxy/T TO O ' x y 



The field enhancement of £ strongly influences the 
field profiles of the heat and charge currents. To dis- 
entangle these effects, we exploit the WF law, which 
states that the ratio n e /Ta is close to the Lorenz factor 
Co = |7r 2 (fcs/e) 2 . In the elemental metals, the WF law 
is nearly universally obeyed at 300 K as well as in the 
impurity-scattering regime below 4 K, while deviations 
are common in between. However, in many interesting 
metals with low carrier densities, n e cannot be measured 
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directly because the charge-neutral term Kb (usually from 
phonons) is comparable in size or larger. 

Recently, a way to separate K e from k using the Righi- 
Leduc effect was introduced. Zhang et al. [1 61 ] have shown 
that the Lorenz ratio may be determined from the ratio 
Ch = K xy /Ta X y. (Essentially, the Righi-Leduc effect 
senses only the electronic entropy current, while filtering 
out the charge-neutral components which do not have a 
Hall response. Since the latter also do not contribute to 
a xy , the ratio of the 2 Hall currents yields the WF ratio. 
The WF-Hall method was tested on high-purity Cu and 
applied to cuprates [II].) CeCoIn 5 is well-suited for this 
method because the 2 Hall conductivities are large. 

The Hall resistivity in CeColn5 was previously re- 
ported [HI], but it is the Hall conductivity that is of in- 
terest here. At each T, we find that the profile of a xy 
vs. H matches that of n xy even when the two curves dis- 
play strong curvature vs. H. The curves of K xy /T and 
C-H&xy are plotted together in Fig. [3] for T < 3 K (Panel 
a), and T > 3 K (b). Let us first note that the curves 
share 2 characteristics rarely seen in Hall experiments. In 
weak H, the curves rise from zero with strong negative 
curvature to produce a knee-like feature. In addition, 
the curvature changes its sign to positive in higher fields; 
both Hall conductivities increase more rapidly than the 
first power in H in strong fields. Further, we note the 
peak anomaly in weak H shown by n xy below 1.5 K. We 
return to these unusual features later. 

At each T, K xy /T and Cncr xy may be scaled together 
over the entire field range by adjusting Ch- We empha- 
size that Ch is an H- independent scaling parameter (oth- 
erwise, it does not make sense to discuss scaling between 
R-xy/T and o xy ). In view of the pronounced nonlinearity, 
the close match between the 2 field profiles is strong evi- 
dence that the WF law is valid with a field-independent 
Lorenz number. The inferred values of Ch are plotted 
in Fig. 01 Between 2 and 10 K, Ch / {ks / e) 2 is close to 
the Sommerfeld value 7r 2 /3, but seems to deviate slightly 
downwards below 2 K. 
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IV. SEPARATION OF ELECTRONIC AND 
NON-ELECTRONIC HEAT CURRENTS 

We next determine K e (T,H) and Kb(T,H) in Eq. [TJ 
Using the values of Ch in Fig. [H we convert the measured 
a into K e (T, H) via 

K e (T,H) = Ta(T,H)C H (T). (2) 

Subtracting the curve of K e (T,H) from k(T,H) at each 
T, we finally determine KbiT, H), which is plotted in Fig. 

Eh.. 

At low T, Kb is found to be strongly H dependent. In 
general, Kb is the sum of the spin-excitation conductivity 
k s and the phonon conductivity K p h, viz. Kb(T,H) = 
k s (T,H) + K p h(T). When spin-disorder scattering of 
phonons is important, an applied field generally leads to 
an increase in K p h because H suppresses spin disorder, 



FIG. 3: Comparison of the thermal Hall conductivity K xy /T 
(solid curve) and the Hall conductivity o xy (open circles) 
scaled by an //-independent Ch at T < 3 K (Panel a) and 
T > 3 K (Panel b). Below 1.5 K (Panel a), K xy /T displays a 
peak anomaly in weak H. Both Hall currents are electron-like 
in sign. 



which is opposite to what is observed. Consequently, we 
identify all the field dependence with the spin-excitation 
term k s (T,H). Figured shows that k s accounts for a 
large fraction of Kb between 5 K and 1 K. At the lowest 
T (0.5-1.5 K), the curve of Kb falls to a floor value at 
the field scale H S (T), which is observable as a break- in- 
slope. In the phase diagram in Fig. QJi, H S (T) is seen 
to lie close to the I/II boundary (solid triangles). With 
the present data, we cannot determine H s above 1.5 K. 
Nonetheless, the trends of the curves in Fig. [T^, suggest 
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FIG. 4: The Lorenz number Lh obtained by scaling K xy /T to 
a xy in Fig. [3] Ch is plotted in units of (fcs/e) 2 . The dashed 
line is the Sommerfeld value 7r 2 /3. 



that, throughout the region I up to 5 K, Kb is greatly re- 
duced from its zero- field values when ff reaches the I/II 
boundary in the phase diagram (Fig. []Jl). 

We compare our results with Ref. [Tj], where the 
Lorenz number was found by estimating the phonon term 
K p h from measurements on a La-doped sample [18]. The 
authors compute the difference k(H,T) - K ph (T), 

assuming K p h to be strictly ff-independent. The in- 
ferred Lorenz number C e ff is reported to be suppressed 
by ~ 30% in a 10-T field [13], in contrast with our ff- 
independent Lorenz number. The field dependence of 
jC e f f in Ref. (l7j | comes from identifying Kb entirely with 
Kph which is assumed ff-independent. However, we note 
that the strong ff dependence of L e fj [ItJ is mainly 
observed below our temperature range. Differences at 
higher T may stem from identifying bosonic contribution 
through experiments on doped samples [l9| . as opposed 
to the Hall response here. In both cases, the Lorenz ratio 
attains the Sommerfeld value above 4 K. 



V. SPIN DEGREES AND CHARGE 
TRANSPORT 

In a conventional magnet, an external H raises the 
magnon dispersion energy which reduces the spin-wave 
population and their thermal conductivity. In the region 
I of CeCoIns, the uniform susceptibility \ is strongly 
enhanced, but conventional long-range magnetic order 
seems to be absent. However, in heavy fermions, spin- 
ordered states involving the local moments in the / bands 
are widely postulated. Incipient spin ordering may exist 
above T c in CeCoIns (see Broholm [2fJ). Although our 
analysis is guided by the known properties of conven- 
tional spin waves, a more exotic kind of spin ordering is 
not precluded, and k s may derive from spin-excitations 
in unconventional spin-ordered states. Because of hy- 
bridization between the / and s-p-d states and large spin- 
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FIG. 5: (a) The curves of Kb — K s + k p h vs. H obtained 
by subtracting K e from the observed k, at each T (Kb is not 
obtained below H C 2 because p — 0). At 0.5, 0.75 and 1 K, the 
H dependence shows a kink at H = H S (T) at the boundary 
between I and II (Fig. QJ.)- The curve of Kb vs. T at 6 T is 
shown in Fig. [lb. (b) Comparison of a xy (bold curves) with 
the quantity a 2 BZ (thin curves). Above the field Hz(T), 
the 2 quantities match over nearly 2 decades with one scaling 
constant Z = 1 x 10 -7 cm 3 /C. Below Hz, however, the scaling 
is spoilt by an "excess Hall current". 



orbit coupling, the spin excitations will strongly scatter 
the charge carriers. We write k s — ^c s v\ where c s is the 
heat capacity of the spin excitations, v the average ve- 
locity and A their mean free path. As the local moments 
align with H, the spin excitation population n s ~ c s de- 
creases steeply in field. We interpret the curves of K b in 
Fig. as the sharp field-suppression of n s at low T. 
While evidence for heat currents carried by spin excita- 
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tions have been reported for low-dimensional oxides [2l| , 
a distinguishing feature in CeCoIn 5 is that changes in 
k s strongly affect the charge currents, which we describe 
next. 

The curves in Fig. [5^ reveal that the charge-neutral 
conductivity Kb decreases strongly with increasing H . 
This trend is opposite to that in the electronic conductiv- 
ity a. As the former decreases, the latter rises in almost 
direct proportion. [The ratio of Kb(H) evaluated at H 
= and 14 T, n b (0) / K b (U) -1.5 and 2.7 at 5 K and 2 
K, respectively. These ratios match the corresponding 
ratios /o(0)/p(14) at the same T in Fig. [5b.] Panel b of 
Fig. [5] shows the steep increase of a 2 and a xy with H . As 
discussed above, the sharp decrease in Kb with H reflects 
a decrease in the density of spin excitations n s . Hence 
the correlated increase in a implies that spin excitations 
are the dominant scattering mechanism of the carriers at 
these temperatures. The steep increase in I in high field 
is caused by field-suppression of the spin excitations. Full 
suppression of this scattering channel, attained when H 
reaches H s , leads to the unusually large I in the region 
II. 



VI. MAGNETORESISTANCE AND CURRENT 
RATIO Tla 

We have also found that the strong enhancement of i 
forges a link between the unusual MR and Hall effect. 
The observed a xy is the sum of contributions a xy from 
each FS sheet i. Assuming that the lifetime Tj on each 
sheet is dominated by spin-disorder scattering, all Tj fol- 
low the same monotonically rising function of field g{B). 
As a result, we have a xy — Bg(B) 2 , while the conductiv- 
ity a ~ g(B): The Hall current grows in direct proportion 
to the square of the longitudinal current multiplied by B. 

To test this assumption, we compare the curves of 
a xy with a 2 B at low T (Fig. [SJa). In high fields, the 
quantity a 2 BZ (thin curves) can be made to match 
&xy (bold curves) by setting the T-independent constant 
Z = 1 x 10 -7 cm 3 /C. The match is excellent for fields 
above a cross-over field Hz{T). As H decreases below 
Hz, however, a xy increasingly exceeds a 2 BZ. The neg- 
ative curvature (knee) feature described earlier now ap- 
pears as a small "excess" Hall current below Hz- 

This high-field scaling is made more apparent if we plot 
the quantity 



U a (T,B)=a xy /a 2 B, 



(3) 



which measures the ratio of the Hall current and the lon- 
gitudinal current squared (see Fig. [6]). We note that 
T^a = Rh[1 + (tan##) 2 ] deviates from the ordinary Hall 
coefficient Rjj when the Hall angle 9h is large. Remark- 
ably, Fig. [Sjshows that, below 1 K, lZ a is just a constant 
equal to Z, even though both a and a xy are increasing 
with strong curvature. Above 1 K, lZ a deviates signifi- 
cantly from Z as the excess Hall current grows, but only 
for fields H < Hz- Above Hz (arrow), we see that lZ a 
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FIG. 6: Curves of the current ratio lZ a (Eq. [3j at selected 
T. Above Hz(T), lZ a attains a plateau value Z which is T- 
independent at our resolution (arrow indicates Hz at 2 K). 
The steep variation in weak H reflects the excess Hall current 
discussed in Fig. [5b. 



again settles down to the value Z. The constancy of TZ a 
is direct evidence that both the anomalous MR and a xy 
reflect the enhancement in I. As seen in the phase dia- 
gram Fig. [Th., Hz (open circles) lies significantly below 
H s . The simple Hall response determines the high-field 
Hall behavior in the regions I and II. Hence the com- 
plicated Hall response in CeCoIn5 arises solely from the 
excess Hall current which is responsible for the weak-field 
"knee" , but is suppressed above Hz ■ 



VII. QUASIPARTICLES IN 
SUPERCONDUCTING STATE 

These findings have interesting implications for the su- 
perconducting state. In heavy-electron systems, pairing 
mediated by the exchange of spin fluctuations has been 
proposed as the likely mechanism for the SC state. How- 
ever, the evidence to date for spin exchange has been 
indirect. Here, we have exploited the unusually large 
strong H dependence of the tensors fty and <Tjj to show 
that spin excitations in fact provide the dominant scat- 
tering channel at low T. Insofar as pairing likely arises 
from the dominant scattering channel, our results provide 
rather direct evidence for spin-mediated pairing. 

Finally, we comment on the extraordinary peak in k vs. 
H that appears below T c (Fig. [7|). The sharp reduction 
of the peak amplitude in H is very similar to the k vs. H 
curves below T c in untwinned YBa 2 Cu 3 7 (YBCO) [lj]. 
The extreme sensitivity to H is in terp reted as caused by 
scattering of nodal qp by vortices |l2l[T3,[l3|. In YBCO, 
the observation of a large anomaly in K xy that peaks at 
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FIG. 7: Expanded view of fcxx Ki VS. H in the vortex state 

of CeCoIns (T < T c ). The sharp peaks in weak H arise from 
the field suppression of the broad peak in k(0, T) below T c 
(see Fig. QJ>). 

finite field provided early key evidence that the peak in 
k arises from enhancement of £ of nodal excitations in 
a ci-wave superconductor. Similarly, the low-field peak 
reported here in K xy at 0.5 and 0.75 K (Fig. [3^) con- 
firms that the cusp-anomaly in k xx is electronic in origin. 
The steep increase in I below T c implies that the nodal 
qp does not experience the intense scattering from spin 
excitations. The close similarity between CeCoIn5 and 
YBCO suggests that a steep enhancement of £ below T c 
associated with nodal quasiparticles may be generic to 
electronic-mediated pairing with d-wave symmetry. 

VIII. DISCUSSION 

We have exploited the unusually large Righi-Leduc ef- 
fect in CeCoIns to determine the Wiedemann- Franz ratio 



of its charge carriers at low T. As shown in Fig. [3l the 2 
Hall conductivities K xy and a xy are strongly non- linear, 
with a change-in-sign of the curvature occuring as H in- 
creases from to 12 T. Remarkably, over a broad inter- 
val of T, the 2 quantities may be scaled together using 
an //-independent Lorenz parameter Ch(T). We find 
that Ch is weakly T dependent and close to the Som- 
merfeld value (it 2 /3)(kB /e) 2 . The strict insensitivity of 
Cr to field allows the electronic heat conductivity K e to 
be determined unambiguously. On subtracting K e from 
the observed k xx , we uncover a large background charge- 
neutral term Kb that is field sensitive. As H increases, Kb 
falls while a rises in proportion. This implies that spin 
excitations are the dominant scatterers of the electrons. 

The transport picture that emerges is that, throughout 
region I in zero H (Fig. QJi), the electrons are strongly 
scattered by spin excitations. Moreover, the spin excita- 
tions contribute the dominant share of the charge neutral 
thermal conductivity Kb, which accounts for ~ 45% of the 
observed k xx at T c . In a finite H||c, the density of spin 
excitations is strongly suppressed. This leads to 2 cor- 
related effects. The neutral heat term Kb is suppressed, 
while the 3 electronic currents cr, a xy and K e grow in 
proportion, as a result of strong enhancement of £. The 
trend in Kb suggests that the full suppression of spin scat- 
tering is attained when H — > H s . An interesting scaling 
relationship between a 2 and o xy is found. (Two recent 
findings related to this work are Refs. [22l[23|.) 
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